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ABSTRACT 

The existence of predominantly cold non-baryonic dark matter is unambiguously demonstrated by several 
observations (e.g., structure formation, big bang nucleosynthesis, gravitational lensing, and rotational curves 
of spiral galaxies). A candidate well motivated by particle physics is a weakly interacting massive particle 
(WIMP). Self-annihilating WIMPs would affect the stellar evolution especially in the early universe. Stars 
powered by self-annihilating WIMP dark matter should possess different properties compared with standard 
stars. While a direct detection of such dark matter powered stars seems very challenging, their cumulative 
emission might leave an imprint in the diffuse metagalactic radiation fields, in particular in the mid-infrared 
part of the electromagnetic spectrum. In this work the possible contributions of dark matter powered stars (dark 
stars; DSs) to the extragalactic background light (EBL) are calculated. It is shown that existing data and limits 
of the EBL intensity can already be used to rule out some DS parameter sets. 

Subject headings: dark ages, reionization, first stars - dark matter - infrared: diffuse background - stars: 
atmospheres 



1. INTRODUCTION 

Several independent observations provide compelling ev- 
idence for an accelerated expanding universe with a matter 
content dominated by a non-luminous ("dark") component 
(ACDM cosmology). Relevant parameters for this so-called 
ACDM cosmology have been measured with high accuracy 



with WMAP (Komatsu et al. 201 1 ) and baryonic acoustic os- 
cillations surveys and will be further refined with upcoming 
Planc k da ta. Large sky surve ys (e.g., SDSS Abazajia n et al | 
|2009| 2dF |Colless et al. 2003 ) and numerical sim ulations like, 
e.g., the Millennium Run ( |Springel et al.|2005| l point toward 
a convincing picture of large-scale structure formation within 
the cosmological concordance model. A promising particle 
candidate for the dark matter (DM) content of the universe 
is the WIMP (weakly interacting massive particle; see e.g., 



Jungman et al. 1996 Bertone et al. 2005 for a review arti- 
cle). Such particles can naturally have annihilation or decay 
channels leading to standard model particles. 

The epoch where the first stars form (z ~ 30) is not yet 
observable with today's astronomical instruments. The cir- 
cumstances and mechanisms of stellar genesis are still top- 
ics of ongoing analysis and rely on sophist icated numerica l 
simulations (for a review see ,e.g., Bromm & Larson 2004). 
Recently, studies discussed the impact of WIMP dark matter 
on the formation of the first stars (|Spolyar et al.|2008 Iocco 



|et al.|2008||focco 2008 ; F reese et al.|2008a| >. During the epoch 
of reionization of the universe, DM may have affected the for- 
mation and evolution of stars. Assuming that self-annihilating 
particles provide the dark matter content of the universe, this 
new source of energy injection into the first stars may alter 
their properties. The energy injection from self-annihilating 
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WIMPs can delay or even prevent the nuclear hydrogen burn- 
ing. 

In principle, two mechanisms could lead to dark matt er ac- 
cretion into a star. The first one, investigated by |Spolyar| 
|et al.| ( |2008[ ), is adiabatic contraction (AC) where additional 
dark matter from outside the forming first star is gravitation- 
ally pulled along with accreting baryonic gas onto the proto- 
star. Due to the very low surface temperature of the Dark 
Star compared to its enormous mass there should be no or 
very little radiative feedback mechanisms preventing further 
accretion. The other possibility to replenish the dark matter 
inside a star is capture of WIMPs via scattering. The effi- 
ciency of this mechanism is depends on the product of the 
elastic scattering cross section between baryonic and dark 



matter and th e surroundi ng dark matter density ( jlocco et al. 
[20081 |loccol|2"00"8l [Freese et al.| [2008a). Both mechanisms 
can lead to generic properties of the DS: low surface temper- 
atures (~ 5 000 - 10000 K), high luminosities (~ 10 6 L ) and 
presumably longer lifetimes than conventional Pop III stars. 
After a "dark phase", the star evolves as a normal zero-age 
main-sequence star. 

Direct detection, even of a Dark Star cluster, is a challeng- 
ing task which may only be possible with future instruments 
like the upcoming J ames Webb Space Telescope (JWST; 



Gardner et al. 
mistic model 

for the search of emission in the early universe is to probe 
the di f fuse m etagalactic radiation field (MRF), see, e.g., [Raue] 
et al.| ( |2009| ). The optical to infrared part of the local MRF 



2006) and under the a ssumption of an opti- 
ackrisson et al.| 2010 >. A different approach 
iuie 



is also known as extragalactic background light (EBL; for a 
review see, e.g., Hauser & Dwek 2001). Its main contribution 
comes from integrated starlight and ther mal dust emissions of 
all cosmic epochs (|Kneiske et al.|[2002 20041 IStecker et al. 
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JFranceschini et al |2008[ JPrimack et al.|2008HGilmore| 
|2009||Finke et al.||2010HKneiske & Dole|2010) . This 



fact makes the EBL a unique probe for the integrated star for- 
mation history of the universe. There are different types of 
observational approaches to measure the EBL. Direct obser- 
vations e.g., with the DIRBE instrument on board the COBE 
satellite suffer from prominent foreground emission, like zo- 
diacal light, that is caused by scattered sunlight b y dust in the 
zodiacal cloud, and diffuse galactic radiation (Hauser et al. 
1998). Lower limits to the EBL are derived from integrated 
galaxy number counts which are availa ble up to a redsh ift * 2 
from Hubble Space Telescope (Madau & Pozzetti 2000 ) and 
the Spitzer instrument (Fazio et al. 2004). A powerful method 
for obtaining upper limits on the EBL density makes use of 
the spectra from very hig h energy (VHE) y-r ay sources, espe- 
cially blazars (see, e.g., |Mazm & Raue 2007 ). Galaxy number 
counts and upper limits from VHE y-ray observations can be 
used to constrain possible DS scenarios. 

In this paper, the so far unknown contribution of DS to the 
EBL is calculated. Constraints for some DS scenarios are de- 
rived as well as a convenient parameterization to calculate the 
maximum EBL density produced by a set of DS parameter 
values. The paper is organized as follows: In Sect. [2] the 
model calculations for the EBL from DSs are described! In 
Sect. [3] the resulting EBL for different sets of DS parameters 
is calculated and compared with recent data and EBL mod- 
els. The multidimensional DS parameter space is constrained 
and a parameterization for the peak EBL contribution of DS 
is derived. Sect. [4] summarizes the obtained results and com- 
pares them with existing direct and indirect approaches to de- 
tect/constrain DS. 

Throughout this paper a flat Friedmann cosmology is 
adopted with O m = 0.3, SIa = 0.7 and a Hubble constant 
of// = 70km s" 1 Mpc" 1 . 

2. THE DS CONTRIBUTION TO THE EXTRAGALACTIC 
BACKGROUND LIGHT (EBL) 

In order to calculate the EBL produced by dark matter burn- 
ing stars a forward evolution model is used based on the calcu- 
lations of Kneiske et al. (2002 ) and Raue et al. (2009) (see also 
|Dweketal.|1998||Salamon & Stecker|1998| ). In the model, the 
specific intensity of the EBL I v (z) is obtained by integrating 
the specific comoving luminosity density s Y (z) over redshift z 



dt 
d? 



Iv(z) = ^ j dz' s y ,(z 

z 

where s v (z) is given by 

-l)Ul\ 

s v (z)= j dz L v (t(z) - t(z'))p*(z') 



(1) 



dt 

dz 7 



(2) 



L Y (t) is the time-dependent specific luminosity, p*(z) the co- 
moving formation rate of DSs at a given redshift z, and 
V = v(l + z)/(l + z') is the redshifted frequency. The integra- 
tion limit Zmax determines the maximum redshift where dark 
matter burning stars begin to form. Cosmological parameters 
enter through 



1 



H Q (\ +z)E(z) 



E{zf = O r (l + z) 4 + + z) 3 + Q k (i + zf + Q A , 



(3) 
(4) 



as described in, e.g., Peebles ( 1993 1. Further details on the 
method and formulae used here can, for example, be found in 
|Kneiskeet ail ( [2002l ). 

DS spectrum with PHOENIX: — The atmospheres and spectra 
of Dark Stars are modeled with the model atmosphere pack- 
age PHOENIX, version 16. A basic description of the code 



can be found in Hauschildt & Baron ( 1999|) and recent appli 
cations to stellar objects can be found in Short & Hauschildt 
(2009) and in Fuhrmeiste r et al.|(|2010[). The cur rent ve rsion 
16 of the PHOENIX package ([Hauschild t & Baron|2010| ) uses 
the ACES equation of state (Barman et al., in preparation) 
to allow for a temperature range from 100K to above 10 6 K 
and includes a multitude of molecular and dust species. The 
models presented here use the ID mode PHOENIX/ID of 
PHOENIX with spherical symmetry. The model atmospheres 
are computed for radiative and convective equilibrium with a 
mixing length to pressure scale height ratio of 2, with the as- 
sumption that no energy is generated in the atmosphere. The 
abundance of H was set to 0.92 by number (mass fraction: 
0.75) and that of He was set to 0.08 by number (mass frac- 
tion: 0.25) for all models, all other elements (including Li) 
have an abundance of zero in these models. Models with ef- 
fective temperatures from 5 00 K to 7 5 00 K with parameters 
(gravities, masses) taken from Spol yar et al.| ( |2009] ) have been 
computed. For all the models a number of different variants 
by, for example, varying the line profiles or the equation of 
state setup to investigate the effects on the models and the syn- 
thetic spectra have been computed. For each setup discussed 
below, the individual models are relaxed to equilibrium before 
computing high-resolution synthetic spectra. 

For a given effective temperature, the choice of line pro- 
files produces the largest variations in the emitted spectra (see 
Fig. [T]). The two sets of models were computed by (a) us- 
ing Stark profiles for the H lines and (b) using van der Waals 
(vdW) broadening for the H lines. At T eW = 5 000 K the vdW 
broadening is the dominant line broadening process whereas 
at T e ff = 7500K Stark broadening is more important, de- 
spite the absence of free electrons from the light metals, see 
the comparison to a spectrum with solar abundances in Fig. [2] 
The hottest model in Fig. [T] clearly shows t he electronic lines 
of H2 in their UV spectra, these are also present in the DS 
spectra with lower effective temperatures but are much weaker 
due to their high excitation energies. In solar abundance mod- 
els these are overwhelmed by the metal lines in the UV and 
would not be detectable. For the model parameters consid- 
ered, the NLTE (Non Local Thermal Equilibrium) effects are 
small and appear to be insignificant. At higher effective tem- 
peratures this will likely be different, however, this parameter 
range is not considered here. 

As can be seen from Fig. [T]no significant hydrogen ionizing 
radiation is emitted for the stellar temperature range consid- 
ered here. 

DS formation rate: — The formation density of DSs n*(z) can 
be linked to the star formation rate p*(z) for the first stars 
(Pop III). |Spolyar et al.j ( [2008] |2009| ) propose that every Pop 
III star that forms inside the center of an undisturbed dark 
mater halo could establish a dark star phase Af D s. This im- 
plies that the mass formation rate of DS can be written as: 
p*(z) = n*(z)xMo^] The comoving star formation rate (SFR) 

1 This assumption holds true if one specific DS model with mass Mds is 
considered. 
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Fig. 1 . — Synthetic spectra for DS models derived with the PHOENIX code. The model parameters are from top to bottom: T e ff = 1 500K, log(g) = 1.0, 
M = 690 M , T cff = 6000K, log(g) = 0.0, M = 371 M Q , 7" cff = 5 500K, log(g) = -0.5, M = 106M o and r cff = 5 000K, log(g) = -0.5, M = 106 M Q . The 
gravitational acceleration g is given in CGS units. 



of the fi rst stars that form in such environments are calcu- 
lated by |Trenti & Sti avelli (2009) with a variety of radiative 
and chemical (metal enrichment) feedback parameters and for 
different numbers of stars forming. They obtain values for the 
Pop III mass formation rate of 10~ 5 - 10~ 3 M Q year -1 Mpc~ 3 
up to z = 10 depending on the exact model parameters con- 
sidered. Pop III star formatio n rates in a similar range are 
also found for lower redshif ts (Schneid er et al.|2006| |Torna^1 



tore et al.|2007 Maio et al. |2010| . There exist two negative 
feedback mechanisms for Pop III star formation: radiative 
feedback from H2 ionizing photons and metal enrichment 
For DSs, due to their cooler temperatures and longer lifetimes 
than Pop III stars, both mentioned feedback mechanisms can 
be suppressed and thus DS formation can, in principle, be en- 
hanced and prolonged. It is also possible that DS can grow 
to larger masses than usual first stars (Mds > Mp op m) which 
results in a higher SFR p*(z). In summary it can be concluded 
that the SFR of DS suffers from large uncertainties. Therefore 
a reasonable fiducial value is derived from the Pop III SFR 
(10~ 5 M G year -1 Mpc~ 3 ) and a wide range limited by extreme 
models (10~ 7 - 10~ 3 M year~' Mpc -3 ) is explored. 

Hence, as a simplification, a constant star formation rate 
over a certain redshift period is assumed which can be ex- 



2 For a detailed description of these processes please refer to |Trenti"&| 

Stiavelli 1 2009 1 and references therein. 

I II f 



pressed as a mass formation rate in units of M year 1 Mpc 

I SFRNorm for 

\0 else 



P*(z) = 



(5) 



where SFRNorm is a normalization factor, varied in the above- 
mentioned range, z m j n indicates the minimal value of redshift 
Z where Dark Star formation can still occur and z max denotes 
the beginning of the Dark Star formation epoch. The ansatz 
for the SFR used in this paper is strongly simplified and more 
elaborate calculations of the SFR are a vailable (see, e.g., |Greif| 
|& Bromm|2006l|Sandick et al.|201 \) . |Raue et al.| l |2003p iave 
investigated the impact of a wide range of SFRs on the re- 
sulting EBL density and found the position and overall height 
of the peak the dominant factor, while the choice of the shape 
only resulted in a weak change in the resulting peak EBL den- 
sity (factor ~2). Given the range of the parameter investigated 
in this work here (e.g., the overall normalization of the SFR 
ranges over several orders of magnitude), this simplified ap- 
proach for the SFR is sufficient. 

As the duration of the DS-forming period in the universe is 
directly linked to the amount of photons that are emitted by 
Dark Stars the influence of z m ; n and z max on the EBL is also 
explored. The contribution of p*(z) for large z is suppressed 
because of the redshift dilution of the photon field which goes 
as (1 + z)~ 3 and so the value of z max is, in the following, set 
to 30. For Zmin (the end of the Dark Star formation epoch) 
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Fig. 2. — Comparison of the Dark Star synthetic spectra for T e g 
T eS = 7 400K, logfe) = 1 .0, M = 1 M Q . 

values between 5 and 15 are considered. These values are 
in good agreement with assumed Pop III formation periods 
( Schneider et al. 200 6l|Tornatore et al.|2007HTrenti & Stiavelli| 
|2009HTrenti et al.|2009HMaio et al.|2010] i. 

DS luminosity: — Independent of the exact mechanism pow- 
ering the DM burning, models predict a stable phase which 
dominates the total radiative output d uring the DS phase 
HSpolyar et aO2008l llocco et al ||2"0W 
the 



7 500K, log(g) = 1.0, M = 690 M Q with a solar abundances model with the parameters 



During this phase 

the lumino sity is nearl y constant (see e.g . Figure 2 in Spolyar 



et al.|20 09 Figure 4 in Iocco et aLp 008 and Figure 1 i n|Yoon 
et al.|20 08 ). Therefore the following ansatz is adopted for L, 



LMD - Kz')) 







for t(z) 
else 



t(z') < Af DS 



(6) 



with Af D s being the duration of this stable phase (also referred 
to as DS lifetime) and being the specific DS luminosity 
according to its synthetic spectrum (cf. Fig. fTh. For Afos 
ffenin) - t(z msx ) the emissivity calculation reduces to 



s v (z)*L° v At DS J" dz'pAz) 



(7) 



leading to a linear scaling of the resulting EBL with Afos- The 
exact length of the DS lifetime is not constrained and depends 
on various factors, e.g., DM type, DS model, DM halo profile, 



etc. (for an extensive discussion see Zackrisson et al. 2010). 
In this work a wide band of possible DS lifetimes is explored, 
ranging from 10 5 to 10 9 years. 

The total bolometric luminosity of a single DS is connected 
to DM particle properties via: 



dvL° oc 



(8) 



( [Spolyar et al.|2008) . 

Given a certain stellar mass formation rate the DS lumi- 
nosity produced per mass has to be specified in the model. 
The mass range of Pop III stars, and, for this reason, also the 
mass range of DS is not very well constrained but is expected 
to be within roughly 10 up to a few hundred solar masses 



(Abel et al. 2002; Schaerer 2002). The model assumptions 
for the DM burning and the star formation lead to a wide 
spread in the DS luminosity per stellar mass. This luminosity- 
to-(stellar)mass ratio (LMR) of published DS model s ( jlocco 
eTaTI [20081 |Spolyar et aT1[2009l |Freese et al"||2Q10) carTBe 
computed and used as input parameters for the EBL calcula- 
tion presented here and range from ~ 10 2 - 10 5 L / M . 

The influence of DM par ticle properties on the LMR will 
be further discussed in Sect. 13.21 



3. RESULTS 
3.1. Constraining the DS parameter space 
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TABLE 1 
Dark Star parameter range 





L Q /M Q 


Ar D s 




SFRNorm 


min 


10 2 


10 5 


15 
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fiducial 


to 3 


10 7 


10 


lO" 5 


max 


10 5 


10 9 


5 


10- 3 



Note. — The maximum value of LMR is close to the Eddington limit and 
therefore acts as a natural upper boundary for the DS parameter space. 



The EBL contribution of DS is calculated using the meth- 
ods and parameters discussed in the previous section. The 
range of the specific parameter values is shown in table [T] 
A fiducial set of "intermediate" parameter values is also dis- 
played which acts as default when a single parameter is var- 
ied. As can be seen from Eqn. [5] in combination with Eqn. 
[T] & Eqn. [2] the resulting EBL density scales linearly with 
the SFR and the LMR. In Fig. [3] the EBL contribution for 
two different DS models is displayed in comparison with a 
strict lower limit for the guar anteed astrophy sical EBL from 
stars and dust in galaxies (Kneiske & Dole 2010). Comparing 
the resulting DS EBL signatures (red dashed and blue dashed 
curves) with their respective input spectra (see Fig. [TJ it can 
be seen that the spectral features are smoothed by the red- 
shift integration. The peak of the EBL is for both models 
located at wavelengths > 2/ mi which diff ers from the value 
(A ~ lyt/m) for Pop III stars (Santos et al. 2002) as expected 
due to lower effective temperatures. The peak of both EBL 
signatures clearly reaches into the detectable region of the in- 
frared background. This shows that the EBL offers the poten- 
tial to constrain DS parameter space. A DS contribution from 
these models would result into a total EBL which is the sum 
of the lower limit and the DS signature (red and blue lines). 
These EBL shapes are already disfavored by upper limits from 
TeV observations (grey line). 

In Fig. |4]the peak EBL contribution for three values of z m i n 
with respect to varying DS lifetimes is displayed. As a conse- 
quence of Eqn. 17] for DS lifetimes smaller than the formation 
period f(z m i n ) - RZmax) the resulting EBL scales linearly with 
increasing Afos- At higher lifetimes than ~ 10 8 years the in- 
tensity of the EBL is increased to a greater amount as well 
as the peak value of the DS signature is shifted toward lower 
wavelengths (Fig. Bl). This is caused by a residual emissiv- 
ity at lower redshiffs z < z m in as the end of DS formation is 
not the end of DS photon emission. Due to the strong dilu- 
tion of the photon number density with redshift (1 + z)~ 3 the 
most recent emission dominates the EBL contribution. If A?ds 
is short enough the end of DS formation is roughly equal to 
the end of DS emitting photons as the amount of DS drops of 
almost instantly. The dashed lines display a linear relation- 
ship between Afos and the maximum EBL flux. Please note 
that DS with lifetimes Atus — 10 10 years, as displayed in Fig. 
H] and Fig. [5] would be still present in today's universe and 
therefore most likely to be detected. 

The model results can be summarized in a formula giving 
the peak EBL contribution at z — from a DS population 
normalized to the fiducial DS parameters. 



(v/ v ) max =2 x 10~ 5 nW nT 2 sr 1 x J 
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Fig. 3. — Two different Dark Star parameter sets (red-dashed: 7Yjs = 7500K, 
M = 690 M Q ; blue-dashed: 7/ DS = 5000 K, M = 106 M Q ). Both models are 
calculated with SFRiMorm = 10~ 3 , Afps = 10 9 years, z m ; n = 5. G rey markers: 
EBL measurements & limits adopted from Mazin & Raue ( 2007 1; grey: upper 
limits from TeV observations (realistic model) from Mazin & Raue (2007). 
Black: EBL lower limit by Kneiske & Dole] (20 10) . The total EBL shape in 
presence of a DS contribution is the sum of the lower limit and the specific 
DS signature (red and blue lines). 
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Fig. 4. — Peak EBL contribution for the 7 500 K DS model as a function of 
DS lifetime, for SFRN orm the fiducial value of 10~ 5 is used. 

mula offers a conservative estimate as the DS lifetime only 
enters linearly which is true for A?ds up to values as large as 
~ 10 8 years. The possible enhancement of the EBL contribu- 
tion due to longer DS lifetimes is not taken into account here, 
but one can estimate it from Fig. |4]or it has to be calculated 
as described in this work. Comparing lower with upper limits 
on the EBL one finds an allowed EBL contribution from DS 
in the range of 5 - 2 5 nW m~ 2 sr~' for w avelength between 
2 and 10/rni (see e.g. |Mazin & Raue|2007| Fig. 15). Adopt- 
ing this range limits on DS parameters can be derived. For 
example: A DS with M DS = 106 M , L DS = 9 x 10 6 L , a 
DS lifetime Afos = 10 years and minimum formation red- 



5 results in a constraint on the DS formation rate 

-3 



shift z„ 

between 5 x 10~ 4 and 3 x 10~ J M year -1 Mpc~ J . In this way 
several DS parameter sets can be used to constrain specific 
DS scenarios. 

3.2. Implications for dark matter properties 

The dominant factor for the DS luminosity due to DM burn- 
ing is the square of the DM density distribution p x which is 



6 



Maurer et al. 




10 10 

Atj, s [years] 

Fig. 5. — Location wavelength of the peak in the EBL SED as a function 
of DS lifetime for the 7 500 K DS model. For Af D s > 10 9 years the values 
converge toward the intrinsic emission maximum indicating that DSs still 
emit light at z = 0. 

determined by the original DM density profile inside primor- 
dial virialized halos, the change of the profile under the influ- 
ence of contracting baryonic material (adiabatic contraction), 
and repopulation of t he inner cusp of the DM halo due to scat- 
tering processes ( |Freese et al.|2008] |2009] >. The characteris- 
tics (temperature, radius, luminosity) of a DS rely on the effi- 
ciency of the DM luminosity compared to other energy gen- 
eration mechanism like nuclear fusion. In general the more 
DM contributes to the total stellar luminosity of the star, the 
cooler and bigger it grows regardless of the mechanism that 
produces the high DM density inside the star. 

As explained above the DS luminosity scales with p 2 . The 
replenishment of the DM density inside a DS is possible either 
by the gravitational pull of the collapsing gas onto the DS or 
by multiple elastic scattering processes between DM and the 
stellar baryons or a combination of both. Th e LMR of the 
specifi c dark matter powered stars presented in |Spolyar et al. 
(2009) tend to be slightly dependent on (crv) ann /m Ar ; higher 
values result in higher DS masses but at the same time also to 
higher DS luminosities. For the early stages of DS formation 
the LMR is nearly proportional t o ^/(o'v) aim /m y. In the case 



of DM capture, as investigated by Iocco et al. (20081), the DS 
luminosity scales linearly with the elastic scattering cross sec- 
tion cr . It is also shown that DM burning due to WIMP cap- 
ture is more efficient in low-mass stars, but under the assump- 
tion that the prod uct (Top x is of the order 10~ 26 GeV cm -1 
|Yoon et aT] ( |2008[ ) have shown that a great mass range of stars 
can undergo a stable DM burning phase. 

It is further worthwhile to note that, while conservative 
models of the Dark Star contribution to the EBL do not 
strongly constrain WIMP properties such as mass or annihi- 
lation cross section, this picture can ch ange when our prio r 
assumptions are relaxed: As outlined in Spo lyar et al.| (|2009 ), 
the Dark Star luminosity-to-mass ratio - and therefore the 
peak EBL contribution (see Eqn. |9| - will roughly scale with 
the square root of the annihilation cross section and WIMP 
mass as LMR oc V(o"v) ann / m where m x is the mass of the 
WIMP. Assuming cross-sections of the order 10~ 23 cm 3 s _1 for 
TeV scale WIMPs, as have recently been invoked to explain 
PAM ELA measurements of the cosmic -ray positron fraction 
( Cirelli & Strumia 2008 ), at least the "max" model presented 
in Table 1 would already be strongly constrained by exist- 



ing EBL data. In comparison to other indirect dark matter 
detection channels, the Dark Star induced EBL component 
is unique in that it is not sensitive to the exact branching ra- 
tios of the annihilation yields into photons or charged particles 
and their resulting spectra, as all annihilation products (with 
the exception of neutrinos, of course) are trapped and ther- 
malized within the Dark Star. In this regard, the calculations 
presented here can also serve as an independent and comple- 
mentary template for indirect dark matter searches in other 
wave bands. 

Implications also exist regarding the elastic scattering cross 
section of the Dark Matter particles. As capture via elastic 
scattering is a dominant channel for rep lenishing the fuel of 
the Dark Star (see e.g. Iocco et aL]|2008 1, once equilibrium 
has been reached the luminosity of such a star scales linearly 
with the WIMP - proton scattering cross-sectior]^] (hydrogen 
of course being the dominant target for an assumed primor- 
dial mixture). There is considerable discussion ongoing with 
regard to experimental searches for elastic WIMP - nucleon 
scattering in underground experiments. 

With respect to Dark Stars, it is interes t ing th at, again 
assuming the basic setup from |Iocco et al.| ( |2008[ ), a spin- 
dependent WIMP-proton cross section of the order lpb for 
a WIMP with a rest mass of tens of GeV would already re- 
sult in Dark Stars with LMRs very close to the "max" value 
of table [TJ While backing up the constraints down to the 
multiple-femtobarn level from th e non-detection of annihi- 
lation neutrinos from the Sun (Hooper et al. 2009), the as 
of yet definitive missing detection of a Dark Star contribu- 
tion to the EBL therefore reinforces the conclusion that spin- 
dependent scattering can reconcile the DAMA and CoGENT 
measurements of an annual modula tion of event rates ( |Bern-| 
|abei et al.|[2010[ |Aalseth et aLl|2011| i with null-results from 
other experiments only for cases in which the WIMP-neutron 
cross-section is several orders of magnitude grea ter than the 
WIMP-proton cross section, as e.g. discussed in Ullio et al. 
(2001 ). It should however be noted that for the case of light 



Dark Matter particles, there exist windows where a surpris- 
ingly high spin-independent scattering cross section of the 
order 10~ 40 cm 2 is comp atible with the DAMA/L IBRA and 
CoGENT measurements ( |Hooper & Kelso 2 01 1) . This may 
again be a case where Dark Matter annihilation might signif- 
icantly effect the life cycle of the first stars, and should be a 
focus of future work. 

4. DISCUSSION & CONCLUSION 

WIMP dark matter can have an impact on the evolution of 
the first stars. In this work, the detectability of DS generated 
signatures in the EBL is investigated. This approach opens 
a new window to search for DM effects: the near-infrared 
(NIR). The EBL contributions from different DS parameter 
sets have been calculated and, for certain sets of parameters, 
the resulting EBL flux can reach into the detectable range of 
the infrared background. A parameterization is presented to 
calculate the peak EBL contribution for a variety of model in- 
put values. The resulting EBL can be used to constrain certain 
DS scenarios by comparing the calculated peak EBL contri- 
bution of DS with existing upper limits in the EBL density. 
The results of this work can be seen as complementary ap- 
proach to put constraints on DS formation scenarios as inves- 

3 In the case of extended adiabatic contraction being the main mechanism 
of DM replenishm ent the DS luminosity scales only with the stellar mass 
jFreese et al.|2010| . 
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tigated by Schleicher et al. (2009 1 where the influence of DS 
to reionization, the y-ray and neutrino background has been 
investigated. 

|Freese et al.| ( |2010| l recently proposed a scenario in which 
Dark Stars can reach enormously high masses (10 7 M Q ) due 
to very efficient and long-lasting accretion processes of bary- 
onic matter onto the Dark Star, so-called super massive Dark 
Stars (SMDS). Under these conditions very long lived and 
luminous SMDS can be produced. As suggested by the au- 
thors nearly all baryonic matter inside the DM halo the SMDS 
forms in is accreted onto the star. This fact translates in a very 
high SFR which can be used as input for Eqn. [9] to put also 
constraints on the SMDS scenarios. 

As a final note, it should be stressed that the ability to con- 
strain DS scenarios via EBL depends crucially on the im- 
provement of EBL limits and measurements. The upcoming 
James Webb Space Telescope will provide new data on deep 
field galaxy counts, as well as refined direct measurements 
of the infrared background. The next generation of Imaging 



Air Cherenkov Telescopes, like the Cherenkov Telescope Ar- 
ray ( The CTA Consortium 2010| l, will radically improve upon 
existing TeV measurements, hopefully providing for a direct 
mea surement of the E BL - induced cutoff in AGN spectra (see 
e.g. Raue & Mazin |2010[ l. Fully self-consistent models of 
Dark Star formation, internal structure and evolution are also 
demanding to finally tackle the question whether WIMP an- 
nihilation played a role in the first luminous objects in the 
universe. 
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